i = 0
fe.results <- list()
for (year in start.year:end.year) {
  i = i + 1
  fe.results[[i]]<- list()
  print(year)
  # Fit FE
  mod <- glm(paste(dv,
                     "~ factor(",
                   aggloc.var,
                   ")",
                     sep = ""),
               data = dta[dta[,t.var]< year,],
               family = "gaussian")
  
  fe.results[[i]]$fit.oos <- predict(mod,
                                       newdata = dta[dta[,t.var] == year,],
                                       type = "response")
  fe.results[[i]]$actual.oos <- as.matrix(dta[dta[,t.var] == year,
                                                      dv])
  print(year)
}
save(fe.results,
     file=paste(modeldir,"/",
                country,
                "_aggfe_",
                v,
                ".RData",
                sep = ""))